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ABSTRACT 



The properties of the asymptotic AdSs space-times representing flat domain walls (DW's) so- 
lutions of the New Massive 3D Gravity with scalar matter are studied. Our analysis is based on 
I st order BPS-like equations involving an appropriate superpotential. The Brown- York boundary 
stress-tensor is used for the calculation of DW's tensions as well as of the CFT2 central charges. 
The holographic renormalization group flows and the phase transitions in specific deformed CFT2 
dual to 3D massive gravity model with quadratic superpotential are discussed. 
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1 Introduction 



The new massive gravity (NMG) represents an appropriate "higher derivatives" generalization of 
3D Eisntein gravity action: 

1 f '1 

Snmg{9^,<t;k,A.) = — dx 3 V=g{eR + —JC-K 2 (-\Va\ 2 + V{a))} 

K = R^R^ - -R 2 , k 2 = 16vrG, e = ±1 (1) 
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which unlike its 4D Einstein version is unitary (i.e. ghost-free) under certain restrictions on the 
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matter potential V(a), on the values of the cosmological constant A = — ^-V(a*) and of the new 
mass parameter m 2 for the both choices of the sign of the i?-term [1]. It is 1-loop UV finite, but 
power-counting non-renormalizable 3 quite as in the case of 4D Einstein gravity [2], [3]. The NMG 
vacuum (a = const) sector contains two propagating (massive) degrees of freedom (the "graviton" 
polarizations) and as a result it admits a variety of physically interesting classical solutions - 
gravitational waves, black holes, etc. [5]. When matter is added the only known exact solutions [6] 
are certain asymptotically dSs geometries that describe Bounce-like evolutions of 3D Universe. 

The problem addressed in the present paper concerns the construction of a family of flat static 
domain walls (DW's) solutions, i.e. a = a(z) and 

ds 2 = dz 2 + e^ z \dx 2 -dt 2 ) (2) 

of the NMG model (1) for polynomial matter potentials V(a). We are looking for DW's interpo- 
lating between two different AdS% vacua (cr\,Af^), parametrized by the solutions of the algebraic 
equations: 



V\a\) = 0, 2A£, (l + = en 2 V(a\) (3) 



The study of such DW's is motivated by their important role in the description of the "holographic" 
renormalization group (RG) flows [7] and of the corresponding phase transitions in two-dimensional 
QFT "dual" to 3D massive gravity (1). The generalization of the superpotential method proposed 
in ref. [6] allows an explicit construction of qualitatively new DW's relating "old" to the "new" 
purely NMG vacua. Assuming further that the AdSs/CFT2 correspondence [8], [9], [10] takes place 
for the extended NMG model (1) as well, we investigate the changes induced by the counter-terms 
K, (and by the sign factor e) on the structure of the corresponding QFT 2 's (3- function, concerning 
the "m 2 -corrections" to the central charges, scaling dimensions and to its free energy. The example 
of the DW's of NMG model with quadratic superpotential and the phase transitions in its dual 
perturbed CFTi (pCFT-z) are studied in some details. An extension of d = 3 NMG's BPS-like I st 
order system and related to it superpotential to the case of d- dimensional New Massive Gravity 
for d > 3 is introduced in Sect. 7. 

2 Superpotential 

Although the NMG action (1) involves up to fourth order derivatives of 3D metrics g^ u , the corre- 
sponding equations for the DW's (2) are of second order: 

a + atp- V'(a) = 



although there exist controversy claims concerning its super-renormalizability [4]. 
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(4) 



due to a particular form of the higher derivatives /C-term. A powerful method for construction of 
analytic non-perturbative solutions of eqs. (4) consists in the introduction of an auxiliary function 
W(a) called superpotential 4 [6], [12] such that 



2eK 2 W 2 {\ 



K 2 V{a) = 2(W') 2 ( y 1 



2em 2 

2„.,a. k 2 W 2 x 



K 2 W 2 ~ 

Aem 2 



(5) 



where W'(a) = & = ^| etc. The statement is that for each given W(a) all the solutions of the 
first order system (5) are solutions of the eqs. (4) as well. For example, the linear superpotential 
W(a) = Ba (B = const) describes a particular double- well matter potential 



(6) 



for em 2 > and 7, m 2 and A given by 



2B 4 k 2 / 5 2 n 2 0/ B 2 

7= — + — ). ml = 8eB 2 (l + — 

\ m, z / \ J77, 



A = m 



The corresponding DW's solutions of eq. (4) 
\/ 2e ,rr|2 



tanh B 2 



Bk 



em" 



( Z - Z 



cosh (b* 
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2e|m 2 | 



(7) 



have as asymptotics at z — > ±00 two very special NMG - vacua with Xbht = — -t- = — 1 [1] placed 
at two degenerate minima 



a 



a(z — > ±00) = ± 



V2em 2 
Bk 



of the potential (6) and representing two AdS% spaces of equal cosmological constant Af^ 
-2em 2 < 0. 



3 Vacua and Domain Walls 

All the constant a solutions of eqs. (5) are determined by the real roots of the following algebraic 
equations:(a) W'{a* a ) = and (b) W 2 (al) = ^pr~, that describe (a part of) the matter potential 
V(a) extrema. Each one of them defines an AdS% space (i.e. one vacua solution of eqs. (5)) 

ds i = dz * + e - 2e ^7 lz (dx 2 - dt 2 ), A = a,b 

4 it represents an appropriate D = 3 NMG adapted version of the Low-Zee superpotential [11] introduced in the 
context of DW's solutions of D — 5 Gauss-Bonnet improuved gravity 
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of cosmological constant A^j = —K 2 W 2 (a A ) as one can see by calculating the values of 3D scalar 
curvature: 

R = -2(p- hp 1 = 8e(W') 2 (\ - - 6k 2 W 2 (8) 

i.e. R vac = — 6k 2 W 2 ((t* a ) = 6A^.j. Hence the variety of admissible vacua of NMG model (1) is 
defined by the values of the extrema 0* of the matter potential V(a) and by the signs of the 
parameters e and m 2 . For example in the case of quadratic superpotential ^2(0) = Ba 2 + D for 
B > and D/Owe find one type (a) vacuum 0* = of cosmological constant = —k 2 D 2 and 
for em 2 > few type (6) vacua given by 



Depending on the range of values of D (i.e. on the shape of potential V(o)) we have: (1) no one 
type (6) vacuum for D > ( 2 ) ^ type (6) vacua for < D < and (3) 

four type (6) vacua {±\a* + ,±\a* + \} for D < -&f?. 

Note that all the type (b) vacua have by construction equal cosmological constants = 
—2em 2 . We consider as an example the DW's one can construct in the region (2) above, charac- 
terized by the three vacua ±|o+| and a* = 0. Then the two DW's solutions of eqs. (5) connecting 
1 + | (or — 1 0+|) with 0* have the following rather implicit form: 

(V) Q+Q - (a 2 + \(a*_) 2 \Y a - ( ( 0+) 2 - a 2 )~ a+ = e ^<-«o) 

£>K.^a_|_CK_ ck_[_k yj 2em^ a_ n yf 2em^ 

e <p- V0 = ^ ^B— ^ )2 _ j\ StB ^ ((7l)2| +u 2^ 8^— (1Q) 

where we have denoted: 

Dk \~ 1 / Dk 



a + = 1 - rr. x ' a ~ = 1 + 



v 7 2em 2 / V v / 2em 2 

Nevertheless one can easily verify that the corresponding asymptotics (at z — > ±oo) of o~(z) 



aiz)^ 00 a* A -a° A e T2AA V^? z , a (oo) = ±\a* + \, a(-oo) = a* a = 

./ -""' 2 - Kff 



A± =-2em 2 , K a eff = -k 2 D 2 , A a = 1 + 



»4 = V>i) (11) 



indeed coincide with the vacuum data (a A , A^j, A^) that determine the boundary conditions 
for the DW solutions in region (2). Observe that the scale factor e v(z ) has different asymptotic 
behaviour depending on the sign of D: in the case of negative values, i.e. for D G (— ^ 2 ^ m2 , O^j and 
e = — 1, m 2 < 0, we find that: 



e*> = -too, e^ z =e- 2 ^^ z ^oo (12) 

while for e = —1, m 2 < and considering positive values within the interval D £ ^0, ^ 2e ™ 2 ^ we 
have that: 

e 2 v^M- oo, ^=°° e 2 v^ -> (13) 
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As it well known the divergences of the scale factor correspond to AdS^ type of boundaries. 
The regions of vanishing scale factor (which are not curvature singularities) represent null Cauchy 
horizons, where the causal description in the Poincare patch terminates. Therefore our DW's (10) 
define particular asymptotically AdS% ((a)AdSs) spaces: the one with D < (see eqs. (12)) has 
two different boundaries and the other one with D > (see eqs. (13)) has one boundary at z — > oo 
and one null horizon at z — > — oo. Let us also mention that the linear superpotential DW's (7) 
(and more general all the DW's relating two type (b) vacua) in the case e = — 1, m 2 < describe 
(a)AdSz spaces of two boundaries, but in this case of equal cosmological constants A^r = — 2em 2 . 



4 Unitarity and BF - conditions 

The Bergshoeff-Hohm-Townsend (BHT) unitarity conditions [1]: 

m 2 (A$ ff - 2em 2 ) > (14) 

together with the Higuchi bound: 

^ ff <M 2 r (A) = -em 2 + ^ ff (15) 

for the massive spin two field ("graviton") are result of the requirement of perturbative (1 - loop) 
unitarity consistency of NMG model (1) for a = const. When massive scalar field is also included 
one have to further impose the Breitenlohner-Freedman (BF) condition [14] which for D = 3 reads: 

A£, < ml(a* A ) = V"{<T* A ) (16) 

or in its stronger form: 

A* ff <ml(a* A )<0 (17) 

It is convenient to parametrize the effective "vacuum masses" m 2 (a A ) = K 2 W A y A (y A — 2) in 
terms of the scaling dimensions = 2 — y A of 2D field <f> a (x,t) "holographically dual" of a(z) [7] 
(A = a,b): 

. 2eW% A k 2 W 2 \ . „. \e{W' h ) 2 

^ = y ^ = ^wX~^)' yb = y{<Tl) = -^w (18) 

Let us consider the case e = —1, m 2 < 0. Then the above unitarity conditions (14), (15), (17) take 
the following simple form 

K 2 W 2 

< ^ < 2, < y a < 2 (19) 

which imposes restrictions on the values of the parameters of NMG model (1). For example, in the 
case of linear superpotential both vacua satisfy all the unitarity conditions only for specific values 
of the parameter B such that: B 2 < \m 2 \. In the particular case of quadratic superpotential of only 
three vacua ±|cr^| and a* = (i.e. in region (2)) all the vacua are unitary and satisfy the weak BF 
condition, when the superpotential parameters are restricted as follows: k 2 D 2 < 2em 2 and B > 0. 
Therefore the corresponding DW's (10) are interpolating between two unitary vacua. Such DW's 
turns out to also have positive tensions tdw > as it shown in Sect. 6. 
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5 Domain Walls Tensions 



In all the "planar" DW's (2) of NMG model (1) the scalar matter is uniformly distributed (i.e. || = 
0) along the whole x-axis and therefore such DW's have infinite energy. As it well known [13], an 
important characteristics of the gravitational properties of such DW's is given by the values of their 
energy densities ejjw = E £ w (equals of their tensions tdw)- In the case of (a)AdSz geometries it 
is given by [15]: 



t dw = lim -L £ v A [ Lx/2 dx?Tjf>?, i,j = 0, 1 (20) 

L x ^oo L x ^ ± J-L x /2 



where A = ± denote the two z — > ±oo limits (8M)a describing (a)AdSs boundaries or/and 
horizons; v± = ±1 and £ M = (0,£ l ) is time-like Killing vector, orthogonal to both (8M)a - surfaces 
and normalized as ClijC = — 1- The Brown- York "boundary" stress-tensor [16] is defined as 
follows: 

Tif = ~ * 5S »* Q v A (21) 
V-7 hi 

where 7^ are the corresponding "boundary/horizon" (8M)a - metrics: 



ltj{x,t)= lim jij(x,t\z), jij(x,t\z) = e^riij, rjij = diag(+,-) 

The main ingredient of the NMG version of the Brown- York formula (20) and (21) is the improved 
NMG action S^\j G = Snmg + S 9 gh with few "boundary" terms S 9 gh added. They represent an 
appropriate generalization of the Gibbons-Hawking boundary action to the case of NMG model (1) 
recently proposed by Hohm and Tonni [17]: 

S gG H = ~ £ v A [ dxdt^ieK - \fK + hijKv) (22) 

where Kij is the extrinsic curvature of 2D "boundary" surface (8M)a\ f^u is the auxiliary Pauli- 
Fierz spin two field [1] whose "on-shell" form 



2 / 1 

Jjxv = —p [Rfiu - -^9^uR ) , m; v = o, 1, 2 



is used in eq. (22); / = j 1 - 1 fij and K = j^Kij. In the case of DW's (2) one can further apply the 
I st order equations (5) in order to derive the following simple "boundary" form of the improved 
action 5 : 



Snmg(DW) = -- £ v A [ dxdt^W(a) (l + 
K A^± J ( 9m )a V 



K 2 W\a) 
2em 2 



(23) 



Then according to the definitions (21 ) and (22) one easily obtains the corresponding explicit form 
of the "boundary" stress-tensor for the NMG model with scalar matter: 



T*{DW) = -\WV> A ) 1 + ^^)75 (24) 



5 unlike the case of the Einstein gravity where all the flat DW's are of BPS type [18], for the 3D NMG DW's 
one need to use the I st order eqs.(5) in order to prove that the remaining terms in the bulk action represent a total 
derivative. 



5 



which allows us to calculate the values of the DW's tensions: 

-nm = - E ± "aWa \1 + ^ J , T^A = W{a\) (25) 

Note that in the m 2 — > oo limit the above formula reproduces the well know results for a flat DW's 
tensions in 3D Einstein gravity obtained by the Israel's thin wall approximation [13]. 



6 Boundary counter-terms and central charges 

Consider NMG model (1) in the limit of small effective cosmological constant: La 3> G = l p i, where 
\kfff\L 2 A = 1. The AdSs/CFT2 correspondence suggests that each of its vacua (a A , A^j, Aa) de- 
termine the main features of certain CFT2 , "living" on the corresponding 2-D boundaries/horizons 
(8M)a of (a)AdS3 space-times (11). As in the case of 3D Einstein Gravity (i.e. the m 2 — > 00 limit 
of (1)), all the 2D data, namely: central charges, scaling dimensions Aa{ct a ) = 2 — i/a and the 
vacuum expectation values 



<A vac \f t j( X+ ,X.)\A vac >=- 2 _ /^g(DW) =T A + T ct, A = T ren {Al x±=x±f (26) 



< A vac I <& a (x+ , X-)\A vac >, etc. of 2D operators and duals of the 3D NMG model fields 7^ 
and a - can be extracted from NMG classical action (23), appropriately renormalized: S'tvmg = 
&NMG + Snmc w here 



act 
o 



The particular form 6 of the boundary counter-terms (27) we have introduced above is a conse- 
quence of the condition that the vacua NMG solutions are conjectured to describe the vacua states 
\A vac >= \<r A ,Afff > of corresponding (UV and IR) CFT2 J s, which by definition must have van- 
ishing dimensions A^ ac = and energy E^ ac =0 (for planar 2D geometries), i.e. 

< T$(x+,x-) > A = = 7# + T% A (28) 

Note that S ct makes NMG action convergent, i.e. we have S r ^^ IG {DW) = 0, thus providing hints 
that such DW's are stable. 

The central charges ca of these CFT2 , s are given by the normalization constants of the stress- 
tensor's 2-point functions 

<f£ ± (x ± )f£ ± (0)> A =^- (29) 

or equivalently by the coefficients of the inhomogeneous part of the stress-tensor's transformation 
laws: 

<S ( it,(x±)>A=~C (30) 



although all our arguments are based on specific DW's solutions of NMG model(l) it is expected that as in the 
Einstein gravity case, this form of the counter-terms is universal, i.e. itcancels the S^g divergences for larger class 
of solutions. 
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under infinitesimal 2D transformations: x± = x± + £±(£±)- According to the Brown- Henneaux's 
observation [19] the 3D counterparts of these 2D conformal symmetries are given by special 3D 
diffeomorphisms that keep invariant the asymptotic form of the AdS% metrics 7 such that 

6 a±± (x + ,x-) = -!f£. (31) 

As a result the corresponding improved Brown- York stress-tensor T[? n (A) (proportional to 7^) 
gets inhomogeneous terms under these transformations: 

klSi A)-&W*(l + !gg)£ (32) 
which allows to calculate the CFT2S central charges in therms of the NMG vacuum data: 

Consider for example the domain wall solution (10) for quadratic superpotential under the restric- 
tion < kD < \/2em 2 and B > , which for e = — 1 and m 2 < interpolates between two 
vacua, i.e. AdS^s of effective cosmological constants A + = —2em 2 and A a = —k 2 D 2 as one can see 
from its asymptotic form (13). Since in this case we have W(a) > 0, the following identification 
kWa = ~tx takes place. As a consequence the corresponding central charges of the two CFT's 
representing these vacua get the familiar form [1] , [17] : 

ca = ^[1 + ^\, L 2 9 r = ^»l 2 pl (34) 

The same central charge formula turns out to be valid in the more general case of DW's for which 
the superpotential W(a) does not change its sign between the two vacua a* A and for the case of 
"non-unitary" vacua with e = 1 and m 2 > as well. It is worthwhile to mention an interesting fact 
that the DW's tensions (25) can be rewritten in terms of the central charges as follows: 

rw(L+ ' ia) = -i(f-|) (35) 

Observe that the condition of positive tensions, i.e. tdw(L + , L a ) > requires |A + | > |A a | which 
is automatically satisfied in the example discussed above. In the case of "unitary" BHT -vacua 
e = — 1 and m 2 < (i.e. for negative c^'s ) this condition is equivalent to the following restriction 
on the central charges: 

fel>|> 1 < 36 > 

i.e. we have c + > c a . It turns out that such "ordering" of the UV and IR central charges determines 
the direction of the RG flow in the dual 2D pCFT2 as we are going to show in the next section. 




7 but are larger then the AdSs isometry group SO(2, 2) 
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7 Comments on holographic RG flows and NMG's extensions 



The off-critical [a)AdS^/pCFT2 version of the holographic principle relates certain static DW's 
solutions of 3D gravity (Einstein or NMG) with scalar matter to the RG flows in specific deformed 
(supersymmetric) CFTi [7]. These non-conformal QFT<is can be realized as an appropriate per- 
turbations of the ultraviolet (UV) CFT2 by marginal or/and relevant operators $ CT that break 2D 
conformal symmetry to the Poincare one : 

S r p e c n FT 2 (v) = Scft 2 + °(L*) J d 2 x^^(x l ) (37) 

The scale-radial duality [20] allows to identify the "running" coupling constant <r(L*) of pCFT2 (37) 
with the scalar field cr(z) and the RG scale L* with the scale factor e^* 2 -* as follows: L* = l v \e~ Lp l' 1 . 
This identification is based on the equivalence of the "radial" evolution equations (5) and the 
Wilson RG equations for the pCFT^. 

da n . , 2eW'(a)/ W 2 (a)K 2 \ , , r 

It is evident that the zeros of the /^-function (38) a* B coincide with the NMG vacuum of type (a) 
{i.e. W'{o* a ) = 0) or of the type (b) (i.e. W±(cr±) = ^pr-)- We also realize that the anomalous 
dimensions A$ of the operator $ CT (x*) at each critical point: 

y(a*) = 2-A^a*) = -^- . (39) 

da <r=a* 

are nothing but the parameters A± b (a*) and y(a* b ) given by eqs. (18), that determine the asymp- 
totic behaviour at z — > ±00 of the matter 3D bulk gravity field a(z). 

As it well known, when the explicit form of the /3(a) - function is given, say by eq. (38), it pro- 
vides the key ingredient that allows to further derive the free energy and certain thermodynamical 
characteristics of 2D classical statistical model related (in its thermodynamic limit) to the quantum 
PCFT2 in discussion 8 . We are interested in the description of the scaling laws, critical exponents 
and the phase structure of particular pCFT2 dual to NMG model (1) with quadratic superpotential 
in the case the range of its parameters B and D belongs to the region (2). Following the standard 
RG methods (see for example [22]) we find that the singular part of the reduced free energy per 
2D volume F s (a) has the following simple form: 

F s (a) « (a 2 )^ ((a* + ) 2 -a 2 )^ (\(a*_) 2 \ + a 2 ) - (40) 

The critical exponents V A = related to the correlation length singularities £4 « (a — a\) y A at 
each critical point (i.e.the NMG's vacua with A = 0,±) are given by: 

4eB 8eB 8eB 

Dn z a + a- a + KV2em 2 a_KV2em 2 

For a particular choice of the "unitary" (for e = — l,m 2 < 0) DW's (10) and of the coupling 
constant a within the range < a < a* + we have that yo < (i.e. the IR CFT2 with ^ a as 



8 indeed we have to consider the euclidean version of NMG such that the corresponding "boundaries/horizons" of 
{a)H'i are fiat euclidean planes or spheres S 2 . 
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irrelevant operator) and < y + < 2 (i.e. the UV CFTi with <J> CT representing now a relevant 
operator). Therefore such DW describes massless RG flow from UV critical point a* + to the IR 
one ctq = 0. An important characteristics of all the massless flows is the so called Zamolodchikov's 
central function : 

which at the critical points a* A takes the values (34). It represents a natural generalization [23] of 
the well known result for m 2 — > oo limit [7], [20]. According to its original 2D definition [21] it is 
intrinsically related to the /3-function: 



3k V da 

of the pCFT2 dual of the NMG model (1). Taking into account the RG equations (38) we realize 
that : 



dC(a) _ 3 fdaV 

dl ~ ~ AGW{a) \ dl) 



(44) 



and therefore when W(a) > is positive (as in our example) the central function is decreasing 
during massless flow we are discussing, i.e. we have c+ > c a . 

Observe that for a > a+ and for a — > oo the correlation length remains finite due to the 
following "resonance" property ^ + ^- + ^- = 0, specific for the quadratic superpotential we are 
studying. Hence this region of the coupling constant space corresponds to the massive phase of 
the pCFT2, which is described "holographically" by singular DW metrics giving rise to (a)AdS3 
space-time with naked singularity, as one can see from the generic form (10) of our DW's solutions. 
We have therefore an example of phase transition from massive to massless phase that occurs at 
the UV critical point a* + . For the description of such phase transition we need two different NMG 
solutions having coinciding boundary conditions {a\, A^, A^) at their common boundary z — > oo. 

We next briefly discuss the possibility to extend our d = 3 superpotential constructions (4) to 
the case of d > 3 NMG models: 

S=±f d d x^-g{R + -L (V iV " J^ZTf 2 ) ~ - 2 {\\^\ 2 + V(a)) } (45) 

As in d = 3 case the static flat DW's solutions of such d-dimensional NMG model are defined by 

d8 2 = dz 2 + e 2 ^T Hj dx i d^, a = a(z), /3 = _ 1 _ , a = (d-l)(3, (46) 

that leads us to a system of second order equations of the type (4), but with different d-dependent 
coefficients. It is then natural to introduce the following generalization of d = 3 NMG superpotential 
and of the I st order system (5) for arbitrary d > 3 : 



2 

2(d-2)m 2 ) y ' 4(d - 2) m 2 
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which in d = 5 case reproduces the Low-Zee superpotential [11] for the Gauss-Bonnet(GB) extended 
5D gravity. It is worthwhile to mention the well known fact (see ref. [2] for example) that for 
conformally flat solutions (i.e. of vanishing d > 3 Weyl tensor as in the case of the DW's (46)) 
the action of the Gauss-Bonnet-Einstein gravity becomes identical to the d > 3 NMG's one (45). 
Therefore the solutions of the eqs.(48) describe the flat static DW's of the both models. The form 
of the eqs.(48) above makes evident that any given superpotential Wd(cr) describe qualitatively 
different matter potentials Vd{a) depending on the values of d = 3,4,5. Hence the properties of 
DW's solutions of corresponding NMG model's, as well as of the /?d(cr)— functions of their pCFTj 
duals, are expected to be rather different depending on the space-time dimensions. We leave the 
problem of the identification of these QFT^ models and of the geometrical NMG's description of 
their phase structure to our forthcoming paper [25]. 

Let us emphasize in conclusion the advantages of the superpotential method in the study of 
the DW's properties of the NMG models (45) as well as of the holographic RG flows in their dual 
pCFTd models. As we have shown on the example of 3D NMG model (1) with quadratic W(a), the 
DW's solutions provide an important information about the phase transitions in its dual 2D model. 
It is important to note however that although we have recognized many of the ingredients of the 
AdS?,/CFT2 correspondence as central charges, scaling dimensions, free energy,etc. the answer to 
the question of whether and under which conditions such correspondence takes place for the NMG 
model (1) remains still open. The complete identification and the description of all the properties 
of the dual pCFT2 in terms of the NMG model's solutions requires better understanding of the 
apparent "unitarity discrepancy" that relates (1-loop) unitary massive 3D gravity to non-unitary 
CFT2 > s of negative central charges in the approximation of small effective cosmological constants. 
Negative central charges are known to appear in different contexts in the (supersymmetric) CFT2S. 

Q 2 

For example, the classical and semi-classical limits % — > of the central charges c q = 1 — 6-^ of the 
so called minimal unitary Virasoro algebra models (as well as of their N = 1 SUSY extensions) 
are big negative numbers [24]. There exist also families of non-unitary 2D models representing 
interesting statistical mechanical problems, as for example the Lee- Yang "edge singularity" CFT 2 
of c = — Tp We remind these facts just to indicate few directions for further investigations that 
might result in the exact identification of the 2D QFT's duals to 3D New Massive gravity models. 
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